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Abstract. Spontaneous edge currents are known to occur in systems of two 
space dimensions in a strong magnetic field. The latter creates chirality and 
determines the direction of the currents. Here we show that an analogous 
effect occurs in a field-free situation when time reversal symmetry is broken 
by the mass term of the Dirac equation in two space dimensions. On a half 
plane, one sees explicitly that the strength of the edge current is proportional 
to the difference between the chemical potentials at the edge and in the bulk, 
so that the effect is analogous to the Hall effect, but with an internal potential, 
fl^ ' The edge conductivity differs from the bulk (Hall) conductivity on the whole 

^ plane. This results from the dependence of the edge conductivity on the choice 

of a selfadjoint extension of the Dirac Hamiltonian. The invariance of the edge 
' conductivity with respect to small perturbations is studied in this example by 

topological techniques. 
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1. Introduction 



When in a two dimensional device without dissipation, an electric field is turned 
on, a current is induced transversally, with density subject to the Ohm-Hall law 
, J = aE. Here a is the 2 x 2-conductivity matrix and an '■— <^2i defines the Hall 

conductivity. For particles described by a Schrodinger operator, a r tiagnetic field - 
perpe ndicular to the plane is needed in addition to obtain gh l)Avron et all 
However, for more general investigations, a time reversal sy mmetry breaking 
\ term in t he H amiltonian might suffice to produce a nonzero cth (|SemenofA Il984l 

■ iHaldanj . Il988,'l . The constant Dirac operator 



. ^ I with fcrmion mass m 7^ yields a very instructive example. Here c is the velocity 

' of light, a :~ (o'i,cr2), where Oi are, for i = 1,2,3, the Pauli matrices, and V is 

^ \ the 2-dimensional gradient. On R^, the operator (Q) features a zero field Hall effect 

l|Frohlich fc Kerleil Il99l|) with an = isgn(m)£^ l|Redlichl. Il98^ . The interpre- 



tation of <7ff at zero temperatur e as the Chern number of a complex line bundle 
ijThouless et all Il982t iKohmotoi Il985l: lAvron fc SeTIel Il985t) fails, but its quan- 
tisati on can be traced back to the geometry of the Lorentz group l|Leitneii 12004 
I2OO5I) . 

In the present paper, we direct our attention to the Dirac operator JQ) on a 
sample with boundary. In this situation spontaneous edge currents may occur, 
wit hout any exterio r electric or magnetic field. We calculate the edge conductivity 
a'^ ^|HalDerinLll982^ for a natural class of self-adjoint extensions of on the half- 
plane. Here ct*^ is an integer (in units of e^ /h) which differs from zero if the boundary 
condition satisfies a certain sign conditi on. It is shown th at a'^ is, in units of 
e^/h, the spectral flow through the gap ijHatsugal Il993albl) . Robustness is then 
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immediate for sufficiently small perturbations of 

In spite of the absence of an exterior field, the edge conductivity can be related to 
the Hall conductivity in the bulk. For Schrodinger operators i n a magnetic field 
equality of bulk (Hall) and edge conductivity has been shown in l|Kellendonk et all 
[2002; Elbau & Graf, 2002). In our system, the relationship is more subtle, since 
the bulk conductivity is half integral, in contrast to the integral edge conductivity. 




Figure 1. Edge currents in the strip. The electric field opens an 
interval A between the respective chemical potentials on the edges. 
The resulting inequilibrium of charge transport along both edges 
yields a nonzero total current in the sample. 



Hall currents can go through the bulk or along the edges, e.g. due to intercepted 
cyclotron orbits l|Schulz-Baldes et all . 1200(1) . For 1, consider the strip [0, c?] x 
M C with edges of similar type at = and xi = d, respectively, and no bulk 
currents. When an electric field of constant strength i?i > is applied, parallel 
to the xi-axis, the force qEi pushes particles away from one edge to the other, 
changing the respective chemical potentials on the left edge {(J-l) and on the right 
edge (/ifl) correspondingly (Figure^]). For /i^ > /i^, there is a net charge transport 
due to the states with energies contained in the interval A := of width 

|A| = If jA(a^i) denotes the corresponding current density, the net edge 

current is 

(2) JI{A) fjAdx,. 

Jo 

It is related to the voltage jAj/g by 

(3) J,^(A)=:a-(A)^. 

Here the pro portionali t y factor a^(A), given in units of q'^/h, defines the edge 
conductivity l|Laughliil Il98lt iHalnerinL Il982l) . (PJ mimicks the Ohm-Hall law 
j2 = a'' El for := (T21. 

For sufficiently large d, the two boundaries decouple, and to calculate J'^(A) = 
J^(A), one only needs to consider a half-plane geometry. Provided /i lies in the 
spectral gap, with lower gap barrier i^crit, A := (iJcrit,M]- 

Depending on whether the boundary is situated on the left (xi = 0) or on the right 
{xi = d) of the sample, the sign in © has to be adjusted, and this is done correctly 
by imposing sgn((T'^(A)) = sgn{a^). |A| = n — Edit > can be interpreted as 
the amount of energy needed to excite a bulk particle of energy i?crit to a state 
at highest possible energy ^. Therefore Q has the shape of the Ohm-Hall law, 
but here the current is proportional to an interior voltage (instead of to an exte- 
riorly applied one as in the Hall effect). In particular, cr'^(A) is again a conductivity. 

Our paper is organised as follows: In the following section, we introduce the 
self-adjoint boundary conditions for the constant Dirac operator on the half 
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plane. Their effect on the spectrum wiU be hivestigated in Section In Sectional 
we derive the corresponding edge conductivity. Section [3 presents a first stability 
result. 

We would like to thank H. Schulz-Baldcs for helpful discussions. 

2. Boundary conditions 

As noticed above, the magnetic field may be zero if a time reversal breaking term 
in the Hamiltonian is present. We investigate the Dirac operator of massive spin 
^ particles (with q ^ e, the electron charge) where this symmetry is broken by the 
mass term. 

D is a symmetric elliptic operator on the domain ^{D) — C^(K+ x R, C^) of 
smooth functions with compact support vanishing in a neighbourhood of xi — 0, 
but it is not essentially self-adjoint. Since D is not bounded below the Friedrichs 
extension is not available for determining a canonical choice of boundary condition. 
Note that even in the Schrodinger/Pauli case, Dirichlet (Friedrichs) and Neumann 
boundary condition are not necessarily the boundar y condition which represents 
the physical system best fsee lAkkermans et aLlHofl^. where chiral boundary condi- 
tions are suggested). Neither Dirichlet nor Neumann nor chiral provide self-adjoint 
boundary conditions for Dirac operators. Therefore, we choose to determine all 
self-adjoint boundary conditions which respect the symmetry of the problem. 

The physical setup is homogeneous w.r.t. X2, and so is D on "D^D). Fourier 
transform in X2 gives a unitary transform 

Jm 

(4) ($(V'))fc2(a;i) := V'fc2(a;i) with 

V'fc,(xi):=^ / e-''=''''ij{xi,X2)dx2. 



An operator is homogene ous w.r.t. xy, if and only if it is decomposable w.r.t. the di- 
rect integral (see, e.g. lR.eed fc Simon. 1978. chapter XIII.16). Of course, we are 
interested only in those self-adjoint extensions D of D which preserve homogeneity. 
We therefore state 

Proposition 1. The X2-homogeneous self-adjoint extensions D of D are given 
exactly by all (measurable) families D{k2) of self- adjoint extensions of D{k2), where 

(5) D{k2) = c^ifcgf^ + cr2hck2 + a-imc^ 
onV{D{k2)) ^C^{R+,C^)- 

Proof. Being a differential operator (with smooth coefficients), D is a closable op- 
erator. By continuity the closure D is homogeneous, and for closed operators we 
have the equivalence between homogeneity and decomposability cited above. The 
fibres D{k2) of Z) are closed, and C^(R+, C^) is clearly an operator core for D{k2). 
This proves the first part. 

The second part is a standard calculation with the Fourier transform. □ 

For determining the self-adjoint extension s of Djky) for fixed fc a we follow the 
von Neumann theory of extensions (see, e.g.. lReed fc SimonLll975l chapter X.l): 

Theorem 1. The self-adjoint extensions of D{k2) are parametrized by Cz R := 
RU {oo}. The extension Z?^(fc2) is given by the domain 

(6) V{D^{k2)) = [{l)^ H\R+) : w{0) = »C^;(0) 
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where Q = oo is understood to mean w(0) — 0, and denotes the L^-Sobolev space 
of order 1. 

Note that, by Sobolev's embedding lemma, i7^-functions on IR+ are continuous, 
so that v{0) makes sense. Physically, ^ says that at xi = 0, no current perpen- 
dicular to the boundary is allowed. Indeed, ji = evi with the velocity operator 
vi := ^[xi, H] — cai acting on C^. Now the matrix element 



(v w) (Ti = vw + wv = 2^{vw) 



vanishes if and only ii w = i(^v for C G M. 

Proof. The bounded parts do not matter for questions of self-adjointness (they do 
change the parametrization) and we choose units with h — \,c — 1 for this proof 
so that we have to deal with T := D{k2) = '^jT'^ only {j — 1, x — xi). 

Since T is first order differential and elliptic, the adjoint is given by the domain 
T){T*) = Vt^^(M+) (i.e. no boundary conditions). According to von Neumann's 
theorem we have to compute the ±i eigenspaces of T* . Because of ellipticity they are 
given by smooth functions, because of uniqueness they are at most one-dimensional. 
We have 

so that V'(a^) = (b)e~^ for some constants a, 6 G C. Reinserting this into the 
eigenvalue equation yields 

(^) 

which is an easily solvable eigenvalue problem in C^. Pj^ := ^(1 ± (jj) are the 
corresponding eigenprojections. To sum up, the ±i eigenspaces of T* are given by 

X± = PfC^e-"". 

Now we have to find all unitaries — > . Since are one-dimensional, 
all unitaries differ only by a complex number z of modulus 1. li k =/= j then 
CTfeCTj = —ajGk by the canonical anti-commutation relations for Pauli matrices. So, 
<JkPj^ = Pf'^k- Therefore, ak maps A'+ to and vice versa, and it is clearly a 
unitary, so that all unitaries are of the form Uz = za^. 

Again, according to von Neumann theory, to each Uz corresponds a self-adjoint 
extension with domain 

(8) I?(T,) = I?(T)e{(l-t/,)V: 

(9) =V{T)®\{l-zau) 



Note that 
so that 



^^V{Tz)^i,{0) = (l-zfTfe)(g) andPr(g) ^0 

(and ip € i?^, of course). In other words, the possible boundary values V'(O) are 
given by the range of i? := (1 - zak)P^ which is a non-othogonal projection. 
Furthermore, 

Pr(l + zak) = i(l - (T,){1 + z<Jk) = 1 - (1 - zak)P+ 
so that the self-adjoint boundary condition can be equivalently described by noting 
(10) P;{1 + zauMO) = ^ V^(O) = (1 - zafc)P+V(0) 
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which we will use in Section |21 

For 3 = 1 and, say, fc = 3, one computes easily = which is 

nonvanishing so that it spans the one-dimensional space of boundary values '0(0) = 
( ^ ) . So we arrived at 

1 + z 

W = V 

1 — z 

which is a fractional linear transformation in z, and as such maps circles to lines 
or circles. Inserting a few values on the circle |2:| = 1 one sees that it is mapped 
indeed to the line iC, C € □ 

Note that, in principle, the parameter C specifying the boundary condition is 
allowed to vary with fc2 without breaking homogeneity. In the following we restrict 
ourselves to constant C, even though the discussion of the spectrum (except for the 
pictures) goes through in the general case as well. 

3. Spectrum 

Note that D^(fc2) depends continuously on k2 so that, by the standard theory of 
direct integrals, the spectrum of is given by 

(11) specZ?^= spec_D^(fc2)- 

fe2eK 

The spectrum of the fibre operator D(^{k2) is determined in the following: 
Theorem 2. The spectrum of D(^(k2) consists of: 



(2) a gap eigenvalue Eg{k2) = — —j^ under the condition 



(1) a continuous part {E: E'^ > Eb{k2)'^}, where Eb = ^ {hck2)^ 
(bulk part) and 

2Chck2 + (1 - C^)™c2 

1 + C 

(12) ftfc2(C^ - 1) > -2mcC. 

Proof. Again we choose the simplified notation from the proof of Theorem ^ and 
write T = Di^{k2). If E is an eigenvalue of T then E'^ is an eigenvalue of 

(13) T'^ = --^ + kl+7n^ 

Ax" 

We begin with the case i?^ < A:| + m? . The only bounded solutions ij} of = 
E^jj have the form 



(14) ^{X) = e'-^kl+ra--E- 

with arbitrary a, 6 G C. Plugging this into the eigenvalue equation T^j = E-tp gives 
the condition 

(15) Qe (l) = E fl) with 



(16) Qe = iy k2 + m'^ — E'^ai + k2(T2 + mas 

in addition to the boundary condition. Note that 

Qe = -{k^ + - E^) + kl + m^ = E'^ 

and tvQE = so that the matrix Qe has spectrum {±i?} and there is always a 
nontrivial solution. For E ^ we define a corresponding (non-orthogonal) eigen- 
projection Pe := ^{1 + ^Qe) (the case i? = is dealt with easily). All candidates 
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for eigensolutions are within the range of Pe- On the other hand, the bound- 
ary condition in the form Hl()|l requires (1 + Z(T3)i/)(0) = 0. A straightforward 
computation with Pauh matrices results in 

(17) A := Pi {1 + za^)PE = 77; (''^ ~ ^fi + — iwa^) where 



(18) V — E — lA/fcl + m'^ — E'^ + zk2i + zm 



(19) w = k2 + zEi- z^kl+m? - E"^ + im. 

The condition for the existence of a nontrivial eigensolution fulfiUing the boundary 
condition is therefore A = 0, since Pe has one-dimensional range onl y. Closer 
inspection shows w = izv so that u = is the only condition to check. (Note that 
the Pauli matrices form a basis of M(2, C).) 



(20) V = Q E + zk2i zm = i\j k\ + m'^ - E'^ 

(21) di{E + zk2t + zm) = and ^{E + zk^i + zm) > 
From this we get 

(oo\ p tc- SI? 2Cfc2 + m(l - C^) 

(22) E = k2^z — mUz = ; — 

1 + C 

and 

(23) < ^zk2i + zm) = k2^z + m^z = '^^^^—^^^^-^^^^ 

1 -f c 

which proves the claim about the gap spectrum. 

In the case i?^ > fc| -f there are always two bounded solutions ■i/j± of T^i/) = 
Eip, having the form 

(24) ^±{x) = e±-Vs^-fe^™^ 

with arbitrary a±,b± € C, so that we have to define two matrices Qe,± and two 
corresponding projections Pe,±- Together with the boundary condition this gives 
the requirement 

= Pf (1 + zas) (pe,+ (lA + Pe 



which has always nontrivial solutions since this is a linear map — > C^. This 
proves the claim about the bulk spectrum. □ 

Remark 1. For the system on R^, D{k2) lives on M, and its spectrum consists of 
{E: E'^ > Eb{k2)'^} only since the solutions for other energies increase exponentially 
either at a; = 00 or a; = —00. This explains the term bulk spectrum because is 
the configuration space of a bulk system. 

For fixed /C2 the bulk spectrum of our 13^ has a gap {—Eh{k2),Eh{k2))- According 
to ijllj, it is A := (-| m|c^, |to|c^). This is the gap we will be interested in. 

Proposition 2. As k2 varies over (—00, 00), the gap eigenvalue Eg{k2) goes through 
the gap (— |m|c^, |to|c^) if and only if m(^ > 0, i.e. when sgnm — sgn^. 

Proof. If = 1 then the gap condition (|12|l requires m( > 0, and Eg{k2) — Qik2C. 
This gives mQ > 0. 

If > 1 then the gap condition requires k2 > kcrit with kcrit ■= Note 
that kcrit is exactly the value of k2 where the line Eg{k2) hits the hyperbola Ei,{k2)- 
Therefore, Eg goes through the gap if and only if kcrit < 0, which is equivalent to 
mC > 0. 
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Figure 2. Spectrum of Dt^{k2) for different m, C- The thick hnes 
are Eg for m, as indicated, the dotted hnes are Eg for —m, 
The dashed hne indicates kcrit (see Proposition 12)). 



If < 1 then the gap condition requires < kcrit- Therefore, Eg goes through 
the gap if and only if kcrit > 0, which is equivalent to mC > again (note that 

— 1 < in the present case, so that the direction of the inequality changes 
again). □ 

4. Edge conductivity on the half plane 

For the constant Dirac operator 111) oyer R , the bulk conductivity is l|RedlichL 
h9S4lT;ndwig et alllT99ilT.eitneTll2nnil2nn,'j^ 

(25) a^ = isgn(m) 

in units of e^//i. To study the corresponding edge conductivity on the half plane, 
for C G let Z)^ = {Z3^(fc2)}fe2gR be the operator family defined by Theorem^] 

Theorem 3. Let A be the gap of the bulk spectrum of Dq. Then for any nonempty 
subinterval A' C A, the edge conductivity cr'^(A') is, in units of 



(26) a^(A') 



sgn(m) if m( > 0, 
otherwise. 



In particular, ct'^(A') = a'^ does not depend on the choice of A' C A. a'^ is the 
spectral flow through E = of D(^. 

Remark 2. The edge conductivity on the half-plane equals the bulk conductivity 
H25|) on in the sense that a'' is the arithmetic mean value of the two possible 
values for a'^. 

Note that interchanging the roles of xi and X2 amounts to rotating the sample 
by 7r/2 and to multiplying C e M by z in the complex plane. If C ^ 0, this yields a 
proportionality factor S R of sign — sgn(C) , and, in terms of of (, the inequality 
in the gap condition of Proposition |2] is reversed. However, this modification leaves 
(T^ unaffected because of the sign convention used in P|). 

Proof. Wc will proceed in two ways. First, let ipk^i^i) be the normalised eigen- 
functions (dj) of Dc_{k2). Eq. Q yields 

dk2 



(27) J^(A) = ec / (T^fe, I a2 \ ^k,)L^iM^) _ - 

Using V2(k2) = h^^dD^(k2)/dk2 and the normalisation condition, we obtain 

1 dEg{k2) _ 2C 



(28) {-iPk, I CT2 Vfe: 



■2/L2(R+) 



ch dk2 + 1 
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from Theorem 121 shows that {tpk2 I j2 \ ipk2)L^{m+) does not depend on k2, so 
that by ©, 

(29) 4^'(A) (X ^ lAj-i / dh 

e JS(fc2)6A 

with proportionahty factor (j2HJ- But r.h.s. of l|29|l is just the absolute value of 
the inverse of the slope of the line Eg{k2). Taking Proposition [21 into account, we 
conclude H26|l . For the last statement, rewrite Q as 

r{A) = eT{lA{H)v2), 

where V2 = -^1x2,0] =: \d2D and 1a{D) being the spectral projection of D onto 
A. T2 is the trace per unit volume in direction X2 for homogeneous operators A, 
defined as 

(30) T2{A) = ^ J^A{k2)dk2, 

where A{k2)dk2 = ^A^~^, and tri is the ordinary trace in direction xi (includ- 
ing the spin-trace over C^). Now approximate IIa/|A| by 5' for a switch function 
g € C°= (l) (denote 1 := R U |±oo j) with g' > 0, suppg' C A, 5(00) = 1, 
g(-oo) = (see, e.g.. lKellendonk et allEool . Then 

J^(A) = er ilAiD)v2) « e|A|r ig'{D)v2) = \A\^t {g'{D)d2D) . 

in 

Denote by '0fc2 ^ normalised eigenvector for Eg{k2), differentiable in fc2- Then 
t^i{g'{D)d2D){k2) - g'{Eg{k2)){^k2\dkMk2)^k,)L-i^+) 

= g'{Eg{k2))^Eg{k2) = ^giE{k2)), 

whose integral is cr'^(A) as given by 126|l . in units of This proof also shows the 
topological nature of the result. 

□ 

The essential point is that instead of varying the subinterval A' C A but using 
the eigenvalue dispersion explicitely, the second approach keeps the calculation 
quite general by introducing a function g which we allow to vary (while now the 
gap interval is fixed). The topological nature of will enable us to show the 
invariance of cr'^ at least under a simple class of perturbations. 

5. Spectral flow and stability 

One of the most remarkable properties of the integer QHE is its stability w.r.t. 
perturbations (disorder). The simplest case is when the perturbation depends on 
xi only: 

Proposition 3. Let W be a bounded self-adjoint operator on i^(M_|_,C^), inducing 
a homogeneous (w.r.t. X2) bounded operator on i^(M_|_ x R, C). // || W^|| < \m\ then 
the system described by -\- W has the same edge conductivity as one described by 

Proof. First note that being bounded, does not change anything regarding the 
boundary conditions and self-adjoint extensions. Since W is independent of X2, the 
direct integral decomposition of + W is D(^{k2) + W, and therefore the Hall 
conductivity is given by the spectral flow as before. 

Through addition of W, the spectrum of -D(fc2) can change by ±||VK|| only. 
Therefore a gap around in the bulk spectrum remains as long as \\W\\ < \m\. 
In the same way, in the -neighbourhood of Eg{k2) there will be a unique 
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eigenvalue of 15^(^2) + 14^ if toC > 0. Since £"^(^2) goes from below — |m| to above 
|m| or vice versa, the unique eigenvalue in the perturbed system will cross in the 
same direction as long as ||W^|| < |m|. Thus the spectral flow is the same. □ 

Note that W is not restricted to be multiplication by a function. Choosing 
W = mi(a;i)CT3 + V(xi) with bounded (smooth, for simplicity) rni^V allows for 
variable mass and electric potential. 

We now turn to the more general case of perturbations which are periodic in X2 ■ 
Since D is not homogeneous w.r.t. X2 any more, we have to replace Fou rier transform 
w.r.t. X2 as in by Floquet-Bloch analysis w.r.t. X2 (see, e.g. .. Reed fc S imon. 1978, 
chapter XIII. 16). Then, for a periodic operator A on i^(R+ x R, C^), its Floquet- 
Bloch transform A{k2) acts on L^(IR_|_ x [0, L], C^) with A:2-quasiperiodic boundary 
conditions on [0, L], and the trace per unit volume is 

(31) T2{A) = ^ [ tTL2[o.L]A{k2)dk2. 

^'^ J[-ti/L,it/L] 

Note that homogeneous operators are in particular periodic, and that for these, 
Definition (ETJ gives the same trace as ()30|l . 




Figure 3. Spectrum in the first Brillouin zone [— tt/L, tt/L]. 
Dashed and dotted lines have the same meaning as in Figure |21 

Reviewing the spectral results from Section |3 in the framework of the Bloch- 
Floquet decomposition leads to the spectrum shown in Figure H Note how m 
this representation (so called reduced zone scheme) the bands and eigenvalues are 
mapped back periodically to the fc2-interval [— tt/L, tt/L]. 

Now, going through the arguments above we see that a"^ is still given by the spec- 
tral flow, even when computed through the Bloch-Floquet decomposition. There- 
fore, Proposition 13 holds mutatis mutandis. 

For physical applications one would like stability under random perturbations 
describing disorder in a crystal. If W is random we cannot apply the Bloch-Floquet 
decomposition any more. I nstead, one could use te chniques from Non-Commutative 
Geometry as was done in teellis sard et al.L Il99 4l) for the quantum Hall-effect. It 
would be interesting to allow randomness in the boundary condition C as well since 
this would describe surface imperfections. We leave this to future work. 
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